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FOREWORD 


This  work  was  performed  by  the  Statistics  Research  Division,  Research 
Triangle  Institute,  Durham,  North  Carolina,  for  the  Electromagnetic  Hazards 
Division,  U.  S.  Naval  Weapons  Development  and  Evaluation  Laboratory,  Dahlgren 


Virginia,  under  Contract  No.  N178-8730. 


INTRODUCTION 


The  energy  required  to  ignite  electro-explosive  devices  (EED's)  of  a 
given  type  will  vary  from  device  to  device.  It  is  not  possible  to  determine 
the  energy  required  to  fire  a  given  device.  The  best  that  can  be  done  is  to 
predict  the  probability  that  a  randomly  selected  device  will  fire  at  a  speci¬ 
fied  energy. 

An  accidental  firing  of  an  EED  by  electromagnetic  radiation  occurs  when 
the  induced  current  exceeds  the  current  required  to  fire  that  particular  EED. 
The  probability  of  an  accidental  firing  may  be  calculated  from  the  distribution 
of  induced  current  and  the  distribution  of  current  required  to  fire  a  randomly 
selected  device.  This  has  been  discussed  by  Richter  cad  Carson  [reference  1]. 
Thus,  any  analysis  of  accidental  firings  of  EED's  will  require  information  on 
the  distribution  of  currents  needed  to  £i.,:n.  EED’s. 

Approximately  8000  Squib  Mark  1  Model  0  devices  were  tested  by  the  Naval 
Ordnance  Laboratory  (NOL)  in  order  to  determine  the  distribution  of  energies 
required  to  ignite  these  squibs.  The  results  of  that  test  and  the  analysis  of 
the  data  have  been  reported  by  Hampton  and  Ayres  [reference  2].  Most  of  the 
data,  particularly  in  the  tails  of  the  distribution,  were  obtained  according 
to  the  Bartlett  plan.  The  Bartlett  plan  used  specifies  che  testing  of  devices 
at  an  energy  level  until  two  reversals  are  observed.  A  reversal  is  a  non-fire 
where  fires  are  predominant  or  vice  versa.  The  NOL  analysis  indicated  that 
the  log-logistic  distribution  was  not  completely  satisfactory,  but  it  provided 
a  better  fit  than  the  log-normal  distribution.  The  NOL  analysis  was  based  on 
maximum  likelihood  estimates  of  the  probability  of  firings  from  the  Bartlett 
plan  data.  In  a  subsequent  study  [reference  3],  NOL  noted  that  a  bias  in  these 
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estimated  firing  probabilities  exists  where  data  are  collected  by  the  Bartlett 
plan.  In  fact,  it  will  be  shown  that  these  estimates  are  biased  high  by 
approximately  a  factor  of  two  at  extremely  low  firing  probabilities.  It  is 
this  lower  tail  of  the  firing  sensitivity  distribution  which  is  of  most 
interest  in  HERO  analyses  since  most  accidental  firings  will  occur  when  the 
current  required  to  fire  an  EED  is  low.  Thus,  it  appeared  desirable  to 
re-anaiyze  the  Squib  MK  1  Mod  0  firing  data  using  an  unbiased  estimator  for 
the  firing  probabilities. 
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Neither  the  log-normal  nor  the  log-logistic  provides  an  adequate  fit 
over  the  entirfd  range  of  the  Squib  Mk  1.  Mod  0  data  fas— stated  in  reference- 
Using  unbiased  estimates  of  firing  probabilities,  both  of  these  dis¬ 
tributions  fit  the  lower  95%  of  the  firing  sensitivity  data  reasonably  well. 
The  log-normal  gives  a  slightly  better  fit  to  these  data.  It  is  illustrated 
that  the  upper  5%  of  the  firing  energies  generally  contribute  a  negligible 
amount  to  the  probability  of  an  accidental  firing.  These  data  indicate  that 
the  log-normal  distribution  provides  a  satisfactory  approximation  when 
calculating  the  probability  of  an  accidental  firing. 

Since  both  the  log-normal  and  log-logistic  fit  the  Squib  Mk  1  Mod  0 
firing  data  about  equally  well,  use  of  the  log-normal  distribution  is 
recommended  in  HERO  computations  where  the  distribution  of  the  induced 
currents  also  must  be  considered  since  this  leads  to  a  simpler  procedure 
for  estimating  accidental  firing  probabilities. 
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DISCUSSION 

Let  r  represent  the  number  of  firings  and  n  the  total  number  of  EED’s 
tested  at  a  particular  energy  level.  If  EED’s  are  tested  at  a  specified 
energy  level  until  r  firings  occur,  an  unbiased  estimate  of  the.  probability 
of  firing,  p,  is  given  by  (r-l)/(n-l),  for  r  >  2,  see  the  Appendix.  For 
the  Squib  Mk  1  Mod  0  data,  r  =  2  with  the  Bartlett  plan  used  by  NOL.  The 

unbiased  estimate  of  the  firing  probability  is  l/(n-l);  whereas  NOL  used  2/n. 

r  „  1 

The  NOL  estimates  in  the  lower  tail  are  high  by  a  factor  of  ( — which 

r-1  n 

is  approximately  a  factor  of  two  for  r  =  2  and  large  n.  Hence,  the  NOL 
analysis  of  the  firing  sensitivity  distribution  is  based  on  biased  estimates 
of  these  probabilities.  The  biased  maximum  likelihood  estimator  does  not 
recognize  that  testing  was  terminated  with  a  firing.  Therefore,  the  observed 
percentage  of  firings  will  be  high.  On  the  other  hand,  if  a  fixed  sample 
size  of  n  EED's  are  tested  at  a  specified  current,  and  r  fire,  then  the 
maximum  likelihood  estimate,  r/n,  provides  an  unbiased  estimate  of  the  firing 
probability.  The  distinction  is:  in  the  first  case  the  number  of  firings,  r, 
is  fixed  arid  the  number  of  EED’s  tested,  n,  is  a  random  variable;  in  the  second 
case,  the  nunber  of  EED’s  tested,  n,  is  fixed  and  the  number  of  fires,  r,  is 
a  random  variable. 

A  description  of  the  experiment  is  given  in  ref.  [2].  The  data  are 
reproduced  in  Table  1  for  convenient  reference.  Unbiased  estimates  of  the 
firing  probabilities  from  the  Bartlett  plan  (negative  binomial)  data  are  pro¬ 
vided  by  l/(n-l),  where  n  is  the  total  number  of  devices  tested  at  a  given 
energy  level.  Sufficient  detail  was  not  given  in  [2]  in  order  to  determine  how 
the  pilot  test  data  v?ere  collected.  It  was  assumed  that  the  sample  size  was 
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fixed  for  the  pilot  test  data.  The  NOL  estimates  were  all  calculated  by, 
r/n,  where  r  is  the  number  of  fires.  This  estimator  is  unbiased  only  for 
the  binomial  (pilot  test)  data  if  the  number  of  teats,  n,  are  fixed  for 
each  energy  level  in  advance  of  the  test.  An  unbiased  estimate  of  the  firing 
probability  at  the  lowest  energy  level  in  Table  1  could  not  be  obtained 
since  it  was  not  known  to  us  whether  or  not  testing  terminated  with  the 
single  firing  that  occurred. 

The  following  notation  will  be  used: 

x^  =  log  energy  at  the  ith  level 

p^  -  unbiased  estimate  of  the  firing 
probability  at  the  ith  energy 
level 

=  standardized  normal  deviate 
corresponding  to  pi# 

If  the  data  fit  a  log-normal  distribution,  then 


y 


(!) 


where  p  is  the  average  firing  log-energy  and  a  is  the  standard  deviation  of 

X  X 

firing  log-energy.  That  is,  y  versus  x  is  a  linear  relationship  with  intercept, 

-p  /a  ,  and  slope,  l/o  .  Estimates  of  p  and  a  can  be  obtained  from  the  fit 
X  X  X  X  X 

of  the  simple  linear  regression  of  y^  versus  x^. 

If  the  data  fit  a  log-logistic  distribution,  then 


(2) 


where  a  and  b  are  constants  which  can  be  estimated  by  simple  linear  regression. 
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TABLE  1.  Firing  information  for  Squib  Mk  1  Mod  0 


Delivered  energy 
(log1Q  millij oules) 

Bartlett  plan 
Fires  Fails 

Pi  lot 
Fires 

tests 

Fails 

Unbiased 

estimate 

.2528 

2486 

2 

.9996 

.2340 

617 

2 

.9984 

.2151 

460 

2 

.9978 

.1964 

819 

2 

.9988 

.1778 

61 

2 

.9839 

.1596 

118 

2 

.9916 

.1408 

37 

2 

.9737 

.1229 

34 

2 

.9714 

.1045 

86 

2 

.9885 

.0864 

14 

2 

.9333 

.0846 

53 

9 

.8548 

.0678 

9 

2 

.9000 

.0508 

5 

2 

.8333 

.0422 

125 

53 

.7022 

.0009 

67 

124 

.  >508 

-  .0158 

2 

3 

.-•  300 

-  .0191 

2 

23 

.0417 

-  .0357 

2 

31 

.0312 

-  .0405 

1 

46 

.0213 

-  .0531 

2 

281 

.0035 

-  .0696 

2 

344 

.0029 

-  .0862 

2 

234 

.0043 

-  .1024 

2 

441 

.0023 

-  .1186 

1 

1270 

. - 
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Ordinary  logit  and  probit  analysis  techniques  are  not  applicable  to 
the  negative  binomial  (Bartlett  plan)  data  obtained  at  NOL  and  presented  in 
ref.  [2].  An  unweighted  least  squares  procedure  was  used  to  fit  the  log- 
normal  and  log-logistic  distributions  to  the  data.  The  way  these  tests  were 
conducted,  the  data  in  the  tails  of  the  distribution  possess  smaller  variances. 
In  a  weighted  least  squares  analysis,  the  ends  of  the  distribution  would 
dominate  the  analysis  and  tend  to  cause  a  straight  line  fit  regardless  of  the 
type  of  distribution.  Hence,  an  unweighted  least  squares  analysis  was  used 
in  order  to  detect  nun-linearity. 


< 
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RESULTS 

The  unbiased  estimates  of  the  firing  probabilities  are  plotted  in  Figure  1. 
The  normal  deviate,  y^,  is  plotted  against  log-energy.  A  straight  line  here 
would  indicate  that  the  data  fit  a  log-normal  distribution.  Log^p^^l-p. )] 
is  plotted  in  Figure  2  against  log-energy.  A  straight  line  here  would  indicate 
that  the  data  fit  a  log-logistic  distribution.  As  can  be  seen  from  Figures  1 
and  2  and  verified  by  statistical  tests,  neither  the  log-normal  nor  log- logistic 
distribution  provides  a  good  fit  over  the  entire  range  of  firing  energy. 

For  the  HERO  problem,  most  accidental  firings  will  occur  with  EED*s  possess¬ 
ing  low  firing  energies.  Hopefully,  for  any  given  set  of  conditions,  the  proba¬ 
bility  of  an  accidental  firing  will  be  low.  That  iss  the  probability  that  the 
firing  energy  is  less  than  the  induced  energy  —  will  generally  be  iow.  A 
satisfactory  estimate  of  firing  probability  can  be  obtained  by  fitting  the  lower 
part  of  the  firing  sensitivity  distribution.  In  fact,  either  a  log-normal  or 
log-logistic  can  be  used  to  fit  all  but  the  upper  5%  of  these  firing  data.  Using 
the  data  with  estimated  firing  probabilities  below  95%  resulted  in  an  extremely 
good  fit  to  the  log-normal  distribution.  The  non-linearity  (quadratic  effect) 
in  the  lower  95%  of  the  distribution  was  significant  at  only  the  20%  level. 

That  is,  if  the  relationship  is  linear  (log-normal  distribution),  the  amount  of 
non-linearity  observed  here  could  be  attributed  to  chance  alone  approximately  80% 
of  the  time  from  similar  experiments.  The  log-logistic  did  not  fit  quite  as 
well,  its  non-linearity  was  significant  at  the  40%  level.  That  is,  60%  of  the 

— ^The  term  "induced  energy"  is  used  here  for  purposes  of  illustration,  since 
the  firing  probabilities  were  estimated  as  a  function  of  energy.  The  energy- 
induced  results  from  the  current  induced  by  electromagnetic  radiation. 
Mathematically,  it  makes  no  difference  whether  energy  or  current  is  used  here 
as  long  as  a  comparable  characteristic  (current  or  energy)  is  used  for  both 
the  firing  distribution  and  the  distribution  induced  by  electromagnetic 
radiation. 
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time  a  worse  fit  would  be  expected  due  to  chance  variation  if  the  distribution 
is  truly  log-logistic. 

Using  the  log-normal  or  log-logistic  distribution  obtained  from  firing 
probabilities  below  95%  extrapolated  over  the  entire  range  of  firing  energy 
will  generally  result  in  a  negligible  over-estimate  of  the  probability  of 
accidental  firings.  The  observed  firing  energies  in  the  upper  tail  are  slightly 
higher  than  those  predicted  by  extending  the  fit  from  the  observed  firing  proba¬ 
bilities  below  95%.  A  few  examples  will  illustrate  that  the  upper  tail  of  the 
firing  sensitivity  distribution  is  generally  unimportant  in  HERO  problems.  Using 
the  Squib  Mk  1  Mod  0  data  below  the  95%  firing  level  resulted  in  an  estimated 
average  firing  log-energy  of  x^  =  .026  with  an  estimated  standard  deviation  of 
s^  =  .039.  Suppose  for  a  given  set  of  conditions  that  the  estimates  of  the  mean 
and  standard  deviation  of  the  induced  log-energy  are  x^,  and  s^.  It  appears  from 
analyses  of  HERO  research  test  data,  Nrval  Weapons  Laboratory,  that  the  log-normal 
distribution  is  a  satisfactory  approximation  to  the  distribution  of  induced  current. 
Since  induced  energy  is  proportional  to  the  induced  current,  induced  energy  would 
be  described  by  a  log-normal  distribution.  Then,  the  difference  between  induced 
log-energy  and  firing  log-energy,  x^  -  x^,  follows  a  normal  distribution  with 

—  —  T” 2  2* 

an  estimated  mean  of  (x^  -  x^)  and  estimated  standard  deviation  of  +  s^  . 

An  accidental  firing  occurs  when  the  induced  log-energy,  x^,  exceeds  the  firing 
log-energy,  x^.  The  approximate  estimated  probability  of  an  accidental  firing  is 

>  x2)  =  P(xx  -  x2  >  0)  . 

The  approximation  results  from  extending  the  log-normal  distribution  from  the 
lower  95%  of  the  firing  sensitivity  distribution  over  the  entire  range  of  firing 
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2  _ 

energy.  Letting  NLp.,0  J  denote  a  normal  distribution  with  mean,  p,  and 
2 

variance,  a  , 

°°  2  2  00 

P(x1  -  *2  —  0)  =  ^  +  a2)]dx  *=  /  N[0, l]dy  ,  (3) 

0  u 

where 

0  -  (atj  -  x2) 

U  - 

For  large  sample  sizes,  this  probability  can  be  approximated  using  tables 
of  the  standardized  normal  distribution. 

The  above  formula  shows  how  firing  probabilities  can  be  estimated  if 
both  the  distribution  of  log-induced  energy  and  the  distribution  of  log-firing 
energy  are  normal.  If  the  log-logistic  distribution  is  used  for  firing  energy, 
then  a  complicated  numerical  integration  procedure  would  have  to  be  performed. 
Since  both  the  log-normal  and  log-logistic  fit  the  Squib  Mk  1  Mod  0  firing 
data  about  equally  well,  use  of  the  log-normal  distribution  is  recommended 
since  this  leads  to  a  simpler  procedure  for  estimating  accidental  firing 
probabilities . 

It  is  of  interest  to  determine  the  contribution  of  the  upper  5%  of  the 
firing  sensitivity  distribution  to  the  probability  of  an  accidental  firing. 
Consider  the  worst  condition  where  all  EED's  in  the  upper  5%  of  the  distribution 
would  fire  at  the  same  energy  obtained  at  the  95%  level  of  the  distribution. 

For  the  Squib,  Mk  1  Mod  0,  it  is  estimated  that  95%  of  the  EED's  would  fire  at 
a  log-energy  below 


x2  +  1.645  s2 


.026  +  1.645 (.039) 


.09  . 
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To  obtain  an  upper  bound  on  the  importance  of  the  upper  5%  of  the  firing 
sensitivity  distribution,  suppose  that  the  remaining  5%  all  fire  at  a  log- 
energy  of  .09.  The  data  show  that  the  upper  5%  of  the  log-firing  energies 
are  higher  than  .09.  Hence,  fewer  accidental  firings  actually  would  occur. 

Under  this  assumption,  the  importance  of  the  upper  tail  is  magnified.  Even 
under  this  assumption,  Table  2  shows  that  the  contribution  of  the  upper  tail 
to  the  probability  of  an  accidental  firing  is  negligible.  Several  values  of 
"x  were  chosen  to  provide  a  wide  range  of  probabilities  which  might  occur 
for  an  accidental  firing.  The  value  of  s^  =  .100  for  the  standard  deviation 
of  induced  log-energy  was  chosen  as  a  somewhat  typical  value  from  HERO  research 
test  data  collected  at  the  Naval  Weapons  Laboratory,  this  value  corresponds 
to  a  coefficient  of  variation,  (standard  deviation  4-  mean/  x  100%,  of  approxi¬ 
mately  25%  for  induced  energy.  The  probability  that  the  induced  log-energy, 
x^,  exceeded  .09,  P(x^  >  .09),  was  then  obtained  from  tables  of  the  standardized 
normal.  The  probability  of  a  firing  log-energy  of  .09  is  .05  under  the 
assumptions.  Hence,  the  probability  that  an  EED  with  a  firing  log-energy  above 
.09  results  in  an  accidental  firing  is  .05  x  PCx^  >  .09).  This  provides  an 
upper  bound  to  the  contribution  to  the  probability  of  a  HE1.0  accident  from  the 
EED's  with  the  highest  5%  of  firing  energies.  An  examination  of  Table  2  indi¬ 
cates  that,  for  HERO-type  probability  calculations,  it  is  not  of  primary 
importance  for  the  distribution  used  to  describe  EED  firing  sensitivities  to 
fit  real  well  in  the  upper  tail. 

It  is  of  interest  to  compare  the  probabilities  of  an  accident  based  on  the 
log-normal  and  the  lpg-logistic  firing  sensitivity  distribution.  It  appears 
that  a  satisfactory  fit  of  the  data  for  HERO  predictions  is  obtained  from  the 
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TABLE  2.  Contribution  of  the  upper  5%  of  firing  currents 
to  the  probability  of  an  accidental  firing. 


Average  induced  energy 

Approximate  prob.  of 

Upper  bound  for 

(log  millijoules) 

accident  .  , 

P(x1  -  x2  >  0)  — 

prob.  of  accident 

*1 

from  EED's  with 
x2  >  .09 

100 

1.20 

X 

10-1 

1.44 

X 

ID'3 

200 

1.74 

X 

10~2 

9.33 

X 

io“5 

300 

1.16 

X 

10-3 

2.40 

X 

10~6 

400 

3.43 

X 

10'5 

2.40 

X 

h~* 

O 

1 

OO 

— !  Calculated  from  (3)  with  the  distribution  of  firing  energy  extrak olated 
from  the  lower  95%  over  the  whole  range  cf  firing  energy. 

2/ 

—  Assumes  that  all  EED's  with  log-firing  energy  above  .09  would  fire  at 
.09. 


rr* 
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lower  95%  of  the  distribution.  The  comparison  is  based  on  extrapolating 
the  log-normal  and  log-logistic  distributions  from  the  lower  95%  of  the 
distributions.  Table  3  shows  a  comparison  of  the  approximate  probabilities 
of  HERO  accidents  based  on  these  two  firing  sensitivity  distributions.  The 
induced  energy  distribution  was  taken  to  be  log-normal  as  before.  The  pro¬ 
babilities  for  the  log-normal  firing  sensitivity  distribution  were  taken 
from  Table  2.  The  probabilities  for  the  log-logistic  firing  sensitivity 
distribution  were  obtained  by  numerical  integration. 

It  is  of  interest  to  note  that  for  these  data,  the  log-normal  distri¬ 
bution  gives  slightly  higher  probabilities.  One  oDjection  that  has  been 
raised  against  using  the  log-normal  is  that  in  past  analyses  by  NOL, 
ref.  [4],  it  has  predicted  higher  firing  currents  in  the  lower  tail  than 
the  log-logistic.  These  data  indicate  just  the  opposite.  This  illustrates 
the  need  for  having  large  sample  sizes  to  estimate  the  lower  tail  of  the 
firing  sensitivity  distribution.  Apparently,  the  variation  in  the  estimate 
of  the  parameters  of  a  distribution  may  have  more  influence  on  probability 
calculations  than  the  form  of  the  distribution  used.  The  close  agreement 
between  the  probabilities  based  on  the  log-normal  and  log-logistic  firing 
sensitivity  distributions  is  encouraging. 
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TABLE  3.  Approximate  probabilities  of  an  accident. 


Average  induced  energy 
(log  millijoules) 

Firing  sensitivity 

distribution 

Log-normal 

Log-logistic 

-1 

_1 

-  .100 

1.20  x  10 

1.11  x  10 

-  .200 

1.74  x  10"2 

1.53  x  10“2 

-  .300 

1.16  x  10-3 

0.96  x  10~3 

-  .400 

3.43  x  10"5 

2.70  x  10"5 

.iLt'/vx<gg»q^g»r|‘gw  wi  a. ; 
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APPENDIX 


The  so  called  negative  binomial  distribution  is  obtained  when  testing 
is  continued  until  a  specified  number,  v,  of  events  are  observed.  The  pro¬ 
bability  that  a  particular  order  of  fires  and  no  fires  results  in  r  fires 

v  n-r 

after  n  tests  is  p  (l~p)  ,  where  p  is  the  probability  that  an  EED  will 

fire  at  a  particular  energy.  Note  that  the  rth  firing  occurs  on  the  nth 
trial.  The  previous  (r-1)  firings  may  occur  in  any  of  the  previous  (n-1) 
trials.  The  number  of  ways  that  (r-1)  firings  can  occur  in  (n-1)  trials  is 
given  by  (n-1) I /(r-1) I (n-r) I  Then,  the  probability  that  n  tests  are  required 
for  r  firings  is  given  by  the  negative  binomial  distribution 


To  show  that  (r-1) /(n-1)  is  an  unbiased  estimator  of  the  firing  proba¬ 
bility  at  a  given  energy,  it  is  necessary  to  show  that  the  expected  value  of 
(r-l)/(n-l)  is  equal  to  p. 


:(£r)=  Jr£r  p(nlr> 


(n-1): 


r  .n-r 
P  (1-P) 


^  (n-2) i  r-1  /n  Nn-r  £  _  0 

Z  <:r-2)T(n-r)T  p  (1'p)  r  >  2 

n=r 


Let  s  =  (r-1)  and  m  =  (n-1);  then 
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E 


“  (m- 1 ) : 

p  JL  T^rn^rp: 

m-s 


p  U-p) 


m-s 


Each  quantity  in  the  summation  is  of 
are  the  probabilities  of  observing  s 
abilities  of  all  possibilities,  m  = 
E  [(r-l)/(n-l) ]  is  p. 

Note  that  this  estimator  exists 


the  same  form  as  (Al).  These  quantities 
firings  in  m  tests.  The  sum  of  the  prcb- 
s,  s+1,  s  +  2,...,  is  one.  Hence, 

only  for  r  >  2. 
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